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Abstract 

In this paper, we perform a systematic study of particle production and neu- 
trino yields for different incident proton energies E p and baselines L, with the 
aim of optimizing the parameters of a neutrino beam for the investigation of 
#i3-driven neutrino oscillations in the Am 2 range allowed by Superkamiokande 
results. We study the neutrino energy spectra in the "relevant" region of the first 
maximum of the oscillation at a given baseline L. We find that to each baseline 
L corresponds an "optimal" proton energy E p which minimizes the required in- 
tegrated proton intensity needed to observe a fixed number of oscillated events. 
In addition, we find that the neutrino event rate in the relevant region scales 
approximately linearly with the proton energy. Hence, baselines L and proton 
energies E p can be adjusted and the performance for neutrino oscillation searches 
will remain approximately unchanged provided that the product of the proton 
energy times the number of protons on target remains constant. We apply these 
ideas to the specific cases of 2.2, 4.4, 20, 50 and 400 GeV protons. We simu- 
late focusing systems that are designed to best capture the secondary pions of 
the "optimal" energy. We compute the expected sensitivities to sin 2 26\^ for the 
various configurations by assuming the existence of new generation accelerators 
able to deliver integrated proton intensities on target times the proton energy of 
the order of 0(5 x 10 23 ) GeV x pot/year. 



1 Introduction 



The firmly established disappearance of muon neutrinos of cosmic ray origin [|I], |2| 
strongly points toward the existence of neutrino oscillations ||. 

1 On leave of absence from INFN Milano. 



The approved first generation long baseline (LBL) experiments — K2K [|J, MI- 
NOS ||, ICARUS and OPERA — will search for a conclusive and unambiguous 
signature of the oscillation mechanism. They will provide the first precise measure- 
ments of the parameters governing the main muon disappearance mechanism. In par- 
ticular, the CERN-NGS beam|8], |J, specifically optimized for tau appearance, will 
allow to directly confirm the hints for neutrino flavor oscillation. 

In addition to the dominant z/ M — > v T oscillation, it is possible that a sub-leading 
transition involving electron-neutrinos occur as well. In the "standard interpretation" 
of the 3-neutrino mixing, the — > v e oscillations at the Am 2 « 2.5 x 10~ 3 eV 2 
indicated by atmospheric neutrinos is driven by the so-called #13 angle. Indeed, given 
the flavor eigenstates v a (a = e, /j,, r) related to the mass eigenstates v\{% = 1,2,3) 
where v a = 11^^, the mixing matrix U is parameterized as: 

C12C13 S12C13 s 13 e~ lS 

U(9i2, #13, 6*23) $) — I — S 12C23 — c 12-5l3S23e* <5 C12C23 — Si2Si^S23e tS C13S23 | (1 

■S12S23 — Ci2Si3C23e J<5 — C12S23 — ■Si2Sl3C23e* <5 C13C23 

with Sij = sin 6ij and c^- = cos 9ij . 

The best sensitivity for this oscillation is expected for ICARUS at the CERN- 
NGS. Limited by the CNGS beam statistics at low energy, this search should allow to 
improve by roughly a factor 5 (see Ref. |J) the CHOOZ|10[ limit on the 813 angle for 



Am 2 « 3 x 10 3 eV 2 . Beyond this program, new methods will be required in order to 
improve significantly the sensitivity. 

At present, the only well established proposal in this direction is the JHF-Kamioka 
project ]TI[. In its first phase, 5 years of operation with the Super-K detector, it aims 



to a factor 20 improvement over the CHOOZ limit. 

In Ref. we have studied an optimization of the CNGS optics that would allow 



to increase the neutrino flux yield at low energy by a factor 5 compared to the baseline 
r-optimization of the CNGS beam. This would yield an improvement in the sensitivity 
by about a factor two, or equivalently an improvement of a factor 10 compared to 
CHOOZ. 



In this paper, FLUKA|T3], [L4j Monte Carlo simulations are employed to perform a 
systematic study of particle production and neutrino yields for different beam energy 
and baselines. Focusing systems adapted to the low energy range are also investigated, 
to obtain realistic estimates of the achievable neutrino rates. 

The whole procedure is assumed to be detector and accelerator independent. Nev- 
ertheless, to give a first estimate of the needed beam intensities, the ICARUS T3000 
detector size has been assumed as a reference. This is also justified by a recent neutrino 



detector comparison |L5|, following which the 2.35 kton LAr ICARUS fiducial mass is 
equivalent to a detector of approximately 20 kton of steel or 50 kton of water. For this 
reason, we concentrate on the intrinsic electron-neutrino background from the beam 
and do not explicitly calculate other sources of backgrounds that are strongly related 
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to detector performances, such as backgrounds from ir° production in neutral current 
events . We assume that in reality the "best" beam would be complemented by the 
"best" detector and that these backgrounds will only introduce a small correction to 
our sensitivity estimates. In addition, we are primarily interested at this stage in the 
comparison among proton drivers. 

We consider so-called "conventional" neutrino superbeams, in which neutrinos are 
produced by the decay of secondary pions obtained in high-energy collisions of protons 
on an appropriate target and followed by a magnetic focusing system. 

In this kind of beams, the neutrino beam spectrum and its flux are essentially 
determined by three parameters that can be optimized appropriately: 



• the primary proton energy E p , 

• the number of protons on target N pot per year, 

• the focusing system, which focus a fraction of the secondary charged pions and 
kaons (positive or both signs depending on the focusing device). 



In order to simplify the problem, we consider three "classes" of proton energies: 



Low energy: protons in a range of a few GeV. We take as reference the CERN- 



SPL proton driver design|16|] with 2.2 GeV kinetic energy and an "upgraded SPL" 
with similar characteristics but with 4.4 GeV protons; 

Medium energy: we take 20 GeV proton energy, similar to the CERN PS machine 
and the 50 GeV of the JHF facility 0; 



High- energy: we take the highest energies, i.e. the 400 GeV protons like in the 
case of the CERN SPS. 



The purpose of this work is to understand the required intensities for the vari- 
ous proton energies, in other words, is it more favorable to employ low-energy, high- 
intensity proton or high-energy, low-intensity proton machines? 



2 Observing the oscillation — choosing L and E 

In order to maximize the probability of an oscillation, we must choose the energy of 
the neutrino E max and the baseline L such that 

L(km) , 0x 7r 
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However, in order to observe the oscillation, we must at least see the maximum preceded 
by a minimum, given by 



The mass difference indicated by Superkamiokande is given by Am 2 « 2.5xl0~ 3 eV 2 
and lies within the range 1 x fCT 3 < Am 2 < 4 x fCT 3 eV 2 Q . 

Since we are considering one single detector at a given location as target, our 
baseline L for a given source is fixed. This implies that the neutrino beam spectrum 
should not be too narrow, e.g. it should be a relatively wide-band beam in order to cope 
with the uncertainty on the Am 2 . Since the knowledge on this parameter will improve 
with time, in particular with the running of K2K, MINOS and CNGS experiments, 
the range of energies will eventually be limited by the desire to observe the minimum 
and the maximum. In the meantime, we consider the full mass range indicated by 
atmospheric neutrino observation. 

We do not wish to consider baselines longer than the current CNGS baseline of 
L = 730 km and use therefore a range of baselines between 100 km and 730 km. 
Table p] shows the values (in MeV) of E m i n and E max defined above as a function for 
the baselines. The energies range from about 50 MeV to 300 MeV for L = 100 km and 
from 300 MeV to 2400 MeV for L = 730 km. Clearly, we wish to optimize for high 
intensity, low energy neutrino superbeams. 

We note that at these low energies, neutrino interactions are clearly identifiable and 
have generally easily reconstructible final-states. This is an advantage for detector- 
related background suppression. In a detector with the granularity like ICARUS, 
neutral pions can be easily suppressed via energy ionization or by direct reconstruction 
of the two (well separated) decay photons. 

3 The approximate scaling with proton energy 

The understanding of the relationship between primary beam energy and neutrino 
production is the first step toward any optimization. To isolate this relationship from 
the other parameters, it is more convenient to work in the "perfect focusing" approx- 
imation, where all (positively charged) mesons produced within a given solid angle are 
supposed to be focused exactly on the detector directionF]. In our case, we set Pt = 
for all secondary positively charged particles produced within 1 rad. 

All the simulations presented in this section have been performed with the same 
target and decay tunnel geometry. The target is in graphite, density 1.8 g/cm 3 , one 

1 ln order to derive the scaling, we consider only neutrino fluxes. However, in Section |6[ when we 
compute our oscillation sensitivities we shall include both neutrinos and antincutrinos components 
obtained from the full simulation of the mesons of the relevant charge. 




(3) 
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404 
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404 


1213 


606 


1617 


809 


600 


485 


243 


970 


485 


1455 


728 


1940 


970 


730 


590 


295 


1180 


590 


1771 


885 


2361 


1180 



Table 1: Neutrino energies E max and E min (see text for definition) corresponding to 
the maximum and minimum of the — > v x oscillation for various baselines and the 
Am 2 range indicated by Superkamiokande. 

meter long, 2 mm radius; the decay tunnel is 150 m long and 3.5 m in radius. Target 
optimization will not change the general scaling, and the effect of the tunnel length 
on muon neutrino production can be accounted for in first approximation assuming 
that all neutrinos are produced by pions, and that they carry the maximum possible 
momentum. 

The pion production rate at the exit of the target for various incident proton energies 
E p as estimated with FLUKA are shown in Figure p] for energies ranging from 2.2 
GeV up to the 400 GeV of the SPS for the CNGS beam. In order to compare pion 
productions at different proton energies, we divide the spectra by the proton energy E p . 

All normalized spectra have similar shape, with the maximum yield at low energies 
(jpy. k, 500 MeV). Far from the endpoint, all spectra scale approximately with the in- 
coming proton energy. Departures from the overall scaling consist in a slightly different 
shape at low E p , and harder spectra at high E p . 

If the detector is far away in the forward direction, the neutrino event rate as a 
function of neutrino energy can be derived from these spectra by considering that: 

• The neutrino carries a momentum that is 0.43 times the parent pion momentum 

• The Lorentz boost gives a factor proportional to E 2 oc E 2 on the solid angle 

• The neutrino cross section grows approximately as E v (not true at the lowest en- 
ergies, where the quasi-elastic grows more rapidly, however, for the optimization 
of the proton energy, this approximation is adequate.) 

For a given E„, these factors apply independently on the primary proton energy. 
Thus, we expect that the energy scaling observed on the pion production is reflected 
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in the neutrino production. This is shown in Figure 0, where we plot the simulated 
muon neutrino event rate scaled with the primary proton energy. 

The super-position of the curves at the lowest energies is impressive, except for 
one at E p = 400 GeV. Of course, by rising E p the energy-integrated event rate raises 
dramatically because the spectra extend to higher energies, but the event rate for a 
given E v <~ 0.15 x E p is simply proportional to E p . 

This approximate scaling implies that we can define a power factor T for neutrino 
production as the product of the proton energy times the beam intensity: 

T = E p x N pot (4) 



Up to now we did not take into account the neutrino oscillation probability. If the 
aim is that of having a neutrino beam centered on the oscillation maximum, E v and 
the baseline L have to be chosen such as 

sin 2 ^1.27Am 2 -^ w 1, (5) 

thus 

E max oc L. (6) 

How to choose the baseline L to maximize the rate of oscillated neutrino events 
per proton incident on the neutrino target? The neutrino flux grows scales like 1/L 2 , 
like the solid angle. But L and E max are proportional, thus dividing bin by bin the 
simulated neutrino spectra with a factor E 2 one gets the shape of the event rate as 
a function of E max . This has been done for the perfect focusing approximation in 
Figure |3] for various proton energies E p . 

It is evident that the most optimal situation is for E max in the range 200-600 MeV 
and this results holds essentially independently of the proton energy, at least as long 
as we discard the very high 400 GeV energy case. We can go beyond the 1/E 2 . ap- 
proximation and in the following we will use the exact expression for the oscillation 
probability and will compute exactly the number of oscillated events by folding the 
expected neutrino spectra. 

Nonetheless, the situation gives us a large freedom in the choice of proton energy, 
provided that the intensity can be re-scaled accordingly, so that the power factor remains 
essentially constant 

T = E p x Npot ~ const. (7) 
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4 The optimal baseline L for each proton energy 



In order to study the scaling with the proton energy E p and the baseline L, we compute 
the number N e of oscillated — > i/ e events: 

POO 

N e oc J pot / dE v a{E v )(j){E v )P{v^ -> z/ e , £„, L, Am 2 , sin 2 20 13 ) (8) 

where 4>{E U ) is the neutrino flux per proton on target. We normalize this number 
to a 2.35 kton (argon) detector, and assume 100% electron detection efficiency. For 
definiteness, we assume unless otherwise noted that Am 2 = 3x 10 -3 eV 2 and sin 2 2#i3 = 
1 x 10- 3 . 

In order to take into account the effect of focusing (realistic focusing is discussed 
in section we focus ideally all particles with the acceptance of 1 rad and apply a 
constant "focusing efficiency" of 20%, i.e. scale down the rates by a factor 5. We will 
see that this assumption is quite realistic and at least conservative for the low energy 
neutrinos. 

Similarly, one can also estimate the "goodness" of the neutrino flux by computing 
the number of muon charged current events N® cc within the energies E m i n (Am 2 = 
1 x 10" 3 eV 2 ) and E max (Am 2 = 4 x 10~ 3 eV 2 ) (see Table 0) of the oscillation: 

rEmax 

Nice* / dE u a(E u )<p(E u ) (9) 

We choose to normalize N® cc to 10 19 pots and per kton of (argon) target. 

With these two parameters, we can easily compare all the possible options. It is fair 
to note that the number of oscillated events N e assumes that the search is not limited 
by background, which in practice is not the case, as we will show in section || However, 
we prefer to separate the issue of background (which will be discussed in section [7J) and 
concentrate for the moment on the question of the needed proton intensity. As we will 
show, the background turns out to be also quite independent of the proton energy and 
hence it does not enter in the optimization of the proton energy. 

The results as a function of the baseline L are summarized in Table |2|. The config- 
uration corresponding to the minimum X pot is shown in bold. 

We observe that for each baseline there is an optimal proton energy E° pUmal , which 
minimizes the required integrated proton intensity X pot to observe a fixed number of 
oscillated events. This is also visible in Figure || where we plot the integrated beam 
intensity needed to obtain 5 oscillated v e events in a 2.35 kton detector for Am 2 = 
3 x 10~ 3 eV 2 and sin 2 (2^13) = 10 -3 for different beam energies and baselines. The 
point of minimum pots corresponds approximately to the maximum event rate. 

This is because in that point, the secondary pion yield energy spectrum per proton 
on target, is best matched to the neutrino oscillation probability. Conversely, for 
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7.7 • 10 24 
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Table 2: Integrated f M CC events per kton and 10 19 p.o.t within the relevant energy- 
interval N® cc and integrated beam intensity T pot , assuming a constant 20% focusing 
efficiency wrt perfect focusing, needed to obtain N e = 5 events (see text) as a function 
of the baseline L for various proton energies. The configuration corresponding to the 
minimum X pot is shown in bold. 



each proton energy there is an optimal baseline L opt , which maximizes the integrated 
neutrino oscillation probability in the neutrino energy region which corresponds to the 
largest weighted pion yield at that proton energy. 

For a 2.2 GeV proton driver, the optimal baseline L opt is approximately L opt rs 
150 km. For a 4.4 GeV proton driver, it is approximately L opt tv 200 km. For 20 GeV, 
we find L opt ps 450 km. For energies above 50 GeV, the optimal baseline is around 
700 km. 

For actual baselines smaller than the optimal baseline, L < L opt , the neutrino 
oscillation maximum occurs at lower energy and the yield of corresponding pions for 
the given proton energy is lower than in the optimal case. Hence, we need a higher 
intensity to compensate for this effect. 

For actual baselines greater than the optimal baseline, L > L opt , the neutrino 
oscillation maximum occurs at higher energy and the yield of corresponding pions for 
the given proton energy is lower than in the optimal case. Indeed, at some point, the 
optimal neutrino oscillation energy corresponds to a pion energy, which is kinematically 
forbidden for the given incident proton energy. Hence, we need again a higher proton 
intensity to compensate for the kinematical suppression. 

At the optimal baselines L opt the power factors are the following: 

w 2.2 x 1.9 x 10 24 = 4.2 x 10 24 GeV x pot, (10) 

^4.4 ~ 4.4 x 8.2 x 10 23 = 3.6 x 10 24 GeV x pot, (11) 

JSo ^ 20 x 1.0 x 10 23 = 2.0 x 10 24 GeV x pot, (12) 

JSo ~ 50 x 5.2 x 10 22 = 2.6 x 10 24 GeV x pot, (13) 
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J^oo ~ 400 x 1.3 x 10 22 = 5.2 x 10 24 GeV x pot. (14) 

This directly confirms the approximate scaling with J 7 , apart at the lowest and the 
highest proton energies. Strictly speaking, a proton energy of 20 GeV appears to be 
the most economical choice in terms of protons, with a proton economy of about a 
factor 2 compared to the 2.2 and 400 GeV cases. 

Summarizing, the required protons on target for 5 oscillated events at Am 2 = 
3 x 10~ 3 eV 2 and sin 2 2# 13 = 10~ 3 in the 2.35 kton mass located at the optimal baseline 
assuming a 20% ideal focusing are k 2 x 10 24 for 2.2 GeV protons, m 10 24 for 4.4 GeV, 
« 10 23 for 20 GeV, « 5 x 10 22 for 50 GeV and « 10 22 for 400 GeV. 

5 Results in real focusing 

The standard focusing system for all neutrino beams up to now is based on magnetic 
horns. However, new solutions should be envisaged for low-energy, intense beams. The 
angular acceptance of the system has to be large (the average transverse momentum of 
reaction products is of the order of 200-300 MeV, comparable with the total momentum 
for the shortest baselines), and the amount of material within the secondary beam cone 
must be as little as possible, both to preserve the flux and to avoid heating/ damage at 
high intensities. 

In principle, the focusing system should cancel the transverse momentum of all 
secondary particles relative to the direction of flight toward the detector and this 
independently of their momentum p. Indeed, this is the definition of the ideal focusing. 

In practice, a horn system can be designed to focalize a given signed momentum, 
i.e. like p. A FODO system focalizes like \p\. A series of coils should focalize like 1, 
i.e. independent of p. Hence, this last solution seems to be attractive in our situation. 

We have found with the help of full simulations and tracking of the particles in 
the focusing system that (1) for the lowest energy configuration, a series of coils does 
indeed provide a quite optimal focusing and (2) for the higher energy configurations, 
where the coil focusing becomes impracticable, the traditional horn focusing can be 
efficiently used. 

5.1 Details on coil focusing 

One possible solution is to exploit the focusing capabilities of magnetic field gradients, 
such as the fringe field at the ends of a solenoid, or the decreasing field far from a current 
loop. The fact that a particle traveling almost parallel to the axis of a solenoid suffers 
a change of its transverse momentum when traversing the fringe region is not intuitive 
but can be understood by realizing that any change of the longitudinal component of 
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an axially symmetric field is associated through Maxwell's divergence equation to a 
radial component of the field. 

The particle motion in axially symmetric fields can be described on the basis of 
the Busch's theorem, which in turn follows from the conservation of canonical angular 
momentum. The derivation and applications of this theorem can be found in textbooks 
(see for instance Jl8|, jl9fl), here we give a summary to better illustrate our focusing 
system. Suppose we have a particle having electric charge q and total energy E, 
moving in a static, non-uniform axially symmetric field, at a distance r from the field 
axis (z— axis), with an angular velocity <p. Busch's theorem states that 

q 

r P<t> + o - ^ = cons t (15) 

where = / -fm<f) is the azimuthal component of the particle momentum, and $ is the 
magnetic flux linked by a circle of radius r centered on the field axis. 

If the field can be considered constant within the circle of radius r, we have 

^ = nr 2 B z , thus rp^ + ~^ r2 B z = const (16) 

If particles are emitted with zero angular velocity by a source located on the axis 
of an axially symmetric magnetic field, both the angular momentum and \l/ are null, 
thus const = 0. 

When exiting to a region of zero field, according to Busch's theorem the particles 
will have zero angular momentum. This is very appealing remembering that in an 
uniform (solenoidal) field, all the transverse momentum is azimuthal. However, this 
is not the case when the field varies, and part of the momentum can also be directed 
along the radius. The radial motion can also be derived exploiting Busch's theorem, 
but the practical use is not straightforward. 

However, if the variation of the magnetic field is small on the scale of the revolution 
time, the theorem of adiabatic invariance states that the magnetic flux encircled by 
the particle trajectory remains a constant of motion. The trajectory radius can be 
assumed constant over a revolution and given by the usual expression 

qB z 

and from adiabatic invariance follows 

B Z K 7r = = const. 

q 2 B z 

Under this condition, all the rotational motion is converted into the longitudinal one, 
and the particle transverse momentum varies as 

# = ir (17) 

-t>2 -DzO 
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where B z0 , the initial magnetic field and transverse momentum. 

Since the period of the motion is 

2tt _ 2nE 

lu b qB z ' 

the condition of adiabatic motion means that AB Z must be small over distances 

2nEv\\ 2np\\ 
Az = - = — — . 

qB z qB z 

For a pion having a longitudinal momentum of 1 GeV/c, in B=l T, Az is about 20 
m. Is is thus difficult to focus efficiently high energy particles through fringe fields, 
while this method can be very effective for low energy ones. It should be stressed 
that this method applies equally well to positive and negative particles, even though 
the resulting advantage is not great since the positive component is dominant for low 
energy beams. 
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3.1 • 10 22 
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0.019 
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0.31 


6.5 ■ 10 22 
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1.5 • 10 22 
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0.009 


2.2 ■ 10 24 


0.027 
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0.38 
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9.8 • 10 21 
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0.005 
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2.5 


8.0 • 10 21 
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0.003 
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0.010 


2.1 ■ 10 24 


0.38 


5.2 ■ 10 22 


3.0 


6.8 • 10 21 


600 


0.002 


9.8 ■ 10 24 


0.006 


3.1 • 10 24 


0.38 


5.2 • 10 22 


3.7 


5.4 • 10 21 
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0.001 


1.6 ■ 10 25 


0.004 


4.9 • 10 24 


0.35 


5.8 ■ 10 22 


3.8 
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Table 3: Same as Table [| with the focusing system included in the simulations. The 
configuration corresponding to the minimum T pot is shown in bold. It is computed in 
order to obtain a number of oscillated events N e = 5 (see text) for the various baselines 
L and proton energies. 



5.2 Results on horn and coil focusing 

We report in Table |3] our results with full detailed simulations of the focusing systems 
for 2.2 GeV 4.4 GeV, 20 GeV and 400 GeV proton energies. 

The coil method has been applied here to the 2.2 and 4.4 GeV proton beams, where 
the produced pions have small energies (only positive mesons have been considered in 
Table ||). The target has been assumed to be a short (30 cm) mercury target, like 
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in the SPL[16j proposal. The non-uniform magnetic field has been obtained with ten 
circular loops, having a radius of 1 m, positioned from to 14 meters from the target, 
carrying decreasing currents to give a field from 20 T to zero. Examples of particle 
orbits and the central magnetic field intensity are shown in Figure f|. The effect on the 
transverse momentum distribution can be appreciated from Figure |^, where positive 
pions emitted from the target within 1 radian have been considered. The decay tunnel 
is 150 m long. 

For higher energy neutrino beams of 20 GeV and 400 GeV, the traditional two-horns 
system has been used. The calculations presented here refer to a first horn to focus 2 
GeV/c particles, followed by a reflector to focus 3 GeV/c. The horn is placed around a 
graphite target, has a length of 4 meters and a current of 300 kA. The reflector starts 
at 6 m from the target, is 4 m long with a current of 150 kA. Examples of particle 
trajectories can be seen in Figure H and the effect on the transverse momentum in the 
case of a 20 GeV primary beam is shown in Figure [7|. The decay tunnel is 350 m long. 

In the total energy range (0-2.5 GeV), the resulting focusing efficiency varies in 
between 0.5 and 0.2. This motivated our choice in Section |], when comparing all 
possible beam/distance options, where a common energy-independent 20% focusing 
efficiency had been assumed. In reality, the focusing profile is not energy independent, 
and the real situation can be better than this, as can be derived from the comparison 
of Tables |2] and 0. This effect is also visible in Figure |8|. 

Summarizing, the required protons on target for 5 oscillated events at Am 2 = 
3 x 10~ 3 eV 2 and sin 2 2^13 = 10~ 3 in the 2.35 kton mass located at the optimal baseline 
assuming optimized realistic focusing are ~ 1 x 10 24 for 2.2 GeV protons, ~ 4 x 10 23 
for 4.4 GeV, « 5 x 10 22 for 20 GeV and w 5 x 10 21 for 400 GeV. 

The cautious reader can wonder why we have apparently been "conservative" in 
assuming a 20% efficiency for the ideal focusing. The point here is that in the ideal 
case we assume a constant efficiency over the whole meson energy range, while in 
the real focusing case one effectively reaches a situation where a part of the meson 
energy range is focalized with an efficiency better than 20% while other parts have 
lower efficiencies. We have of course optimized focusing for the energy relevant to 
the oscillation. It would however be incorrect to assume that this higher efficiency is 
constant over the full energy range. 

6 A superbeam to Gran Sasso (BNGSQ) ? 
6.1 Finding the location for the source 

As a working hypothesis, we take for granted the existence of LNGS as an underground 
laboratory that can host large neutrino detectors. In particular, we have assumed for 

2 BNGS stands for "Better Neutrinos to GS". 
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Table 4: List of existing laboratories with interesting baselines to the LNGS Gran 
Sasso laboratory. Angle is the space angle relative to the orientation of the LNGS 
Halls. Inclination is the incoming neutrino angle relative to the horizontal plane at 
LNGS. 

this study the existence of ICARUS with 2.35 kton fiducial mass. 

If we take the location of the detector as fixed, the baseline L is determined by 
the location of the neutrino source. We have investigated various potential locations 
within Italy where large ENEA0 or INFN infrastructures are already existing, as shown 
in Table |j. In these locations, the required conditions to host a high intensity proton 
driver could be met and the machine would find other applications in addition to 
neutrino physics. 

In Table (|, the column called "angle" describes the space angle relative to the 
orientation of the LNGS Halls. Inclination is the incoming neutrino angle relative to 
the horizontal plane at LNGS. We note that due to the fortunate orientation of the 
LNGS Halls according to the geographical axis of the Italian peninsula, it is possible 
to find various laboratories within the orientation of the LNGS Halls. 

Indeed, ENEA Aquilone, ENEA Brasimone, ENEA Trisaia, INFN Pavia and ENEA 
Ispra appear with an angle less than 15° with respect to the LNGS Hall direction. This 
is an advantage for the acceptance of higher energy events given the natural longitudinal 
orientation of the detectors. On the other hand, in the case of ENEA Casaccia situated 
near Rome, the angle is 63°. However, for the shortest baseline, we expect the relevant 
neutrinos to have an energy similar to that of most atmospheric neutrinos, so that we 
can argue that the acceptance will not be a problem given the isotropical nature of a 
detector like ICARUS. 



3 ENEA stands for "Ente per le Nuove Tecnologie, l'Energia e l'Ambiente" . 
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Table 5: Neutrino beam parameters for a 270 km baseline experiment. All quantities 
are calculated for < E u < 2.5 GeV 

6.2 The medium baselines (120 < L < 470 km) 

We study the cases ENEA Casaccia (L = 120 km), ENEA Brasimone (L = 270 km) 
and INFN Pavia (L = 465 km). At these distances, the neutrino energy range relevant 
for — > v e search is 0.1-1.0 GeV. We assume 2.2, 4.4 and 20 GeV proton energies. 
Expected neutrino fluxes and rates, obtained assuming ideal (p.f.) and real (coil or 
horn) focusing systems, are reported in Table [|. The figures are normalized to the 
baseline L = 270 km and for < E v < 2.5 GeV. Other baselines can be rescaled 
accordingly. With real focusing, we find about 100(10) v^ip^) CC events/kton per 
10 23 pots for 2.2 GeV protons, about 340(40) v^) CC events/kton for 4.4 GeV 
protons and about 26000(250) v^ip^ CC events/kton for 20 GeV protons. Intrinsic 
v e {p e ) beam contaminations are in the range from 1% to 2% with respect to v^iv^. 

The 2.2 and 4.4 GeV proton energies are the closest to the optimal energies and only 
these two cases are considered in the following. We will consider 20 GeV in section |6l| 
in the context of the CERN-LNGS baseline. 

To appreciate the matching between the neutrino beams and the oscillation prob- 



ability, we show in Figures |9] and ^ the charged current event rate at 120 km as a 
function of the neutrino energy. The dotted lines correspond to the oscillation proba- 
bility (arb. norm.) for a Am 2 = 3 x 10 -3 eV 2 . 

Clearly, these neutrino beams offer optimal condition to study z/^ disappearance. 
Indeed, the maximum of the oscillation is very well covered by the neutrino beam 
and hence it is quite obvious that a very precise determination of the main oscillation 
parameters will be accomplished. We here do not consider this any further. 

We, however, concentrate instead on the u e appearance measurement. Since we are 
in the presence of intrinsic v e background from the beam at the level of l%-2% of the 

component, we can improve our sensitivity by studying the energy spectrum of the 
v e charged current events. This method is more sensitive than simple event counting. 
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In order to estimate the sensitivity, we adopt our standard fitting procedure of the 
various reconstructed event classes (See Ref. [0). We assume that the neutrino and 
antineutrino interactions cannot be distinguished on an event-by-event basis, and hence 
add the v e and v e contributions from the beam. Similarly, the oscillated spectrum is 
calculated by summing both —>■ u e and — > u e oscillations, assuming the same 
oscillation probability for neutrino and antineutrinos. 

In the present study, we considered only the energy distribution of electron events 
and computed the \ 2 as a function of the sin 2 20 13 mixing angle, scanning in Am 2 . 
The 90% C.L. sensivity region is defined by the condition \ 2 > Xmin + 4-6, defined by 
the condition that the actually observed events in the experiment coincide with the 
expected background. In an actual experiment, a simultaneous fit of the muon disap- 
pearance and electron appearance spectra will constrain the Am 2 , sin 2 23 parameters 
and in case of negative result will limit the sin 2 20i 3 within the allowed Am 2 region. 

The results of the Am 2 scans are shown in Figure |ll| for 2.2 GeV and Figure [l^ 



for 4.4 GeV proton energy for the three assumed baselines. The curves correspond to 
5 years running with 2 x 10 23 pots/year equivalent to a continous proton current of 
1 mA and a fiducial mass of 2.35 kton. The assumed protons on target is compatible 
with an accelerator with performances similar to those of the planned CERN-SPL||16||. 



The expected sensitivities represent a great improvement to the CHOOZ limit JTO 
which gives for Am 2 = 2.5 x 1CT 3 eV 2 



(sin 2 28 13 )chooz < 0.14 or 13 < 11° (18) 

Indeed, we find for 2.2 GeV protons: 

(sin 2 26 13 ) BNG s, 120 k m < 0.006 or 13 < 2.2° (19) 

(sin 2 26 13 ) BNGSt270km < 0.015 or 13 < 3.4° (20) 

(sin 2 2e 13 ) BNGS ^ km < 0.03 or 13 < 5° (21) 

and for 4.4 GeV protons: 

(sin 2 29 13 ) BNGS ,i Wkm < 0.0035 or 13 < 1.7° (22) 

(sin 2 26 13 ) BNG s,270k m < 0.006 or 13 < 2.2° (23) 

(sin 2 26 13 ) BNGSAe5k m < 0.02 or 6 13 < 4° (24) 



For comparison, the JHF proposal with the OAB beam gives similar results fTT[ 

(sin 2 20i 3 ) jhf,oab < 0.006 or 13 < 2.2°. (25) 

6.3 The CERN-GS baseline 

The baseline between CERN and GS is 730 km. At this distance, the neutrino energy 
range relevant for — > u e search is 0.3-2.5 GeV. 
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The present CNGS design|| is optimized for v T appearance, thus for a relatively 
high-energy neutrino beam. The 400 GeV/c SPS beam will deliver 4.5 10 19 protons per 
year on a graphite target, made of spaced rods to reduce the re-interaction rate. The 
two magnetic horns (horn and reflector) are tuned to focus 35 and 50 GeV/c mesons, 
with an acceptance of the order of 30 mrad. The decay tunnel length is 1 km. With 
the standard CNGS parameters, the low-energy neutrino flux is low, as can be seen 
from the entries flagged by > in Table |6|. 

In Ref. ||12|| , we have studied a L.E. optimization of the 400 GeV protons of the 
CNGS in order to improve the sensitivity to 6\ 3 . This yielded an improvement of a 
factor 5 in flux at low energy compared to the r optimization. 

Here, we study the 20 GeV proton energy (we call this the PS++) and compare it 
to 400 GeV. Expected neutrino fluxes and rates obtained with real focusing systems 
are reported in Table []. 

The results of the Am 2 — sin 2 2# 13 sensitivity scans are shown in Figure O. The 



curves correspond to 5 years running with 2 x 10 21 or 2 x 10 22 pots/year and a fiducial 
mass of 2.35 kton. 

We find for Am 2 = 2.5 x 10~ 3 eV 2 : 

(sin 2 26 13 )p S++i2xW 2ipot/year < 0.016 or 6»i3 < 3.6° (26) 
(sin 2 29 13 ) PS++!2xW 22 pot/year < 0.005 or 9 13 < 2° (27) 



Clearly, a strongly intensity-upgraded CERN PS booster would provide very inter- 
esting possibilities for the oscillation searches over the CERN-GS baseline. 
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Table 6: Neutrino beam parameters for the CNGS baseline, with E u < 2.5 GeV. The 
' cases correspond to the present CNGS design for target, acceptance and focusing 
system. 
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7 The intrinsic v e background 



As well known, the intrinsic electron (anti)neutrinos in the beam are produced either 
in the decay of muons coming from kaons or pions via the chain 



7T + /K + -> /!+ + , yU + -> z/ e + e + + f M , 



or directly in the three-body i^ e 3 kaon decays. 

When looking for — > z/ e oscillations, this contamination will eventually be the 
limiting factor. It is therefore essential to understand its level and it is also worth 
understanding if the beam design can be optimized to minimize this background. 

Rather than the v e J ratio, we decide to consider the ratio v^j yJUl in order to 
better estimate the effect of the backrgound on the — > v e oscillation sensitivity. 

Which source of v e is relevant to our study? 

Due to the large difference between the ir and \x decay lengths, the electron neutrino 
background depends on the length of the decay tunnel /. This dependence should be 
more evident for low beam energies. 

Kaon production strongly depends on the proton energy, as shown in Figure [T5|, 



where a threshold effect is clearly visible. For the lowest proton energies, 2.2 and 
4.4 GeV, the fraction of kaon relative to pions is a few per mille. Above 20 GeV, it is 
on the order of a little less than 10%. 

At low v e energies, the production is shared by kaon and muon decays, while kaon 
decays alone are responsible for the high energy tail. For low proton beam energies, 
kaon production is much lower and practically all the intrinsic electron contamination 
comes from muon decays. 

A first order estimate of the effect of the decay tunnel length I on the muon-induced 
background can be derived assuming forward decay at each step, and counting for each 
neutrino energy the fraction of parent particles that decay within a path /. The fraction 
D w of 7r decayed after a length I is simply given by 

D w = (l-e-^) (28) 
Muons have to be generated by a pion first, thus 

v / t-v 



= j A^e • (1 - e V J dy 

= 1 - ( - K*~*) ( 29 ) 



The decay lengths A^, A M depend on the energy of the meson. We can fix a neutrino 
energy E u , the same for and v e . For production, we can then assume p n = 
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Ev/0.43. For the u e produced in the decay of a muon, there is no fixed relation, but 
we can take the average of the v e energy in the muon rest frame, getting p^ w E u /0.Q. 
The grand-parent pion had therefore a momentum equal to p' n ~ p^ ~ E u /0.6. 

With these approximations, £)„- and Z) M can be expressed as a function oll/E v . In 
these units, one 7r decay length corresponds to l/E v — 130 m/GeV. 

It is a priori obvious that short decay tunnels reduce the relative probability of 
muon production and decay However, the pion decay yield is also affected. 

To study the background from muons, we consider the statistically significant ratio 
7r decays/ a/ fidecays as the correct estimator for v^j^fvl as a function of the decay 
tunnel length. This ratio is shown in Figure [HJ where the assumption of a fixed v e 
energy in /x decay has been relieved: the hatched band corresponds to 20%-80% of the 
maximum v e energy in the \x rest frame. 

We find that the ratio does not show dramatic variations between and 4 pion 
decay lengths. We therefore conclude that not much is to be gained by reducing the 
length of the decay tunnel. 

We have verified these results directly by the full simulation of the neutrino beams 
for various decay tunnel lengths. Answers are reported in Table [7] for various proton 
energies and decay tunnel lengths. The 6th column shows the expected v e contamina- 
tion relative to the and the last column lists the statistically relevant ratio v^j -^v~ e . 
For the 2.2 GeV proton energy, the ratio v^jv^ varies from 0.3% for I = 20 m up to 
1.7% for I = 150 m, but this happens at a high cost of genuine v^s. The statistically 
relevant ratio v^j ' yfv~ e varies from 0.67 for I = 20 m down to 0.47 for I = 150 m. This 
is a modest loss. We also stress that the naive y/~N scaling is not adequate for an 
appearance experiment where we are looking for few events and hence we conclude 
that genuine rate is more important than a slightly better v e j ratio. 

Remains the issue of the proton energy. Naively, one would expect that the higher 
the proton energy, the higher is the background. We find however that the intrinsic 
electron neutrino background does not strongly depend on the proton energy. 

This was verified directly for various proton energies and baselines. Results of the 
calculation are shown in Table |8|, all normalized to 10 19 pots. We observe that (1) for 
the shortest baselines, the ratio v^/v^ is increasing dramatically with proton energy. 
Accordingly, (2) the ratio v^js/vl decreases. However, we must rescale this ratio to 
take into account the approximate scaling of the number of events with the proton 
energy. Since we expect 

Vtl E n r—— 

" OC OC y/W p (30) 



\l E p 



we consider the rescaled ratio v^j ^Jv e -E v in the last columns of the Table |[ These 



are also plotted in Figure [17] as a function of the baseline L. Numerically, we find that 



the rescaled ratios v^j w v e ■ E p are almost the identical at the optimal baselines of each 
proton energy, so not much is too be gained by varying the proton energy. 



18 



Summarizing, we find that the electron neutrino background in the relevant region 
is not dependent on the proton energy and only determined by the decay tunnel length, 
but that its optimisation is very limited. For the maximum neutrino muon flux, it is 
at the level of the 1% for any of the considered setups. 
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Table 7: Electron neutrino intrinsic background within E v < 2.5 GeV for various proton 
energies and beam optics configurations. For ease of comparison, rates are normalized 
to a baseline L=730 km. 



8 Summary and Conclusions 

In this document, we have performed a two-dimensional scan, varying the beam energy 
and baseline parameters to optimize the conditions for the investigation of #13 driven 
neutrino oscillations in the whole Superkamiokande allowed Am 2 range. We find that: 

• The optimal baselines for 9 is searches are in the range 100-700 km for proton 
energies varying from 2.2 to 400 GeV. 

• The needed beam intensity scales approximately with the inverse of the beam 
energy. 

• In terms of proton economics, the optimum beam energy is around 20 GeV, but 
lower beam energies are appealing for the shortest baselines. 

Realistic focusing system for low and medium baselines have also been studied. In 
this case: 
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Table 8: Relationships between muon neutrino CC events and electron neutrino CC 
events, for 10 19 pots. Event spectra have been integrated over the energy range of 
interest at each baseline L. 

• Focusing efficiencies of 30-50% can be achieved in the energy range of interest. 

The whole procedure is detector and accelerator independent. Nevertheless, to give 
a first estimate of the needed beam intensities, the ICARUS detector has been assumed 
as a reference. 

We can draw the following observations: 

• a 2.2 GeV or 4.4 GeV high-intensity proton machine (i.e. a la CERN-SPL) is well 
matched to a baseline in the range 100-300 km. It is not matched to a baseline 
of 730 km (i.e. CERN-LNGS). 

• a 20 GeV machine is best matched to a baseline of 730 km (i.e. CERN-LNGS). 
However, an integrated intensity in the range of 10 23 pots are required, which 
is about two orders of magnitude higher than the intensity deliverable by the 
current CERN-PS in a reasonable amount of time. 

• a 400 GeV energy is reasonably matched to a baseline of 730 km (i.e. CERN- 
LNGS). For the 400 GeV, the required intensity is in the range of 10 22 pots, 
which is about one order of magnitude higher than the intensity deliverable by 
the current CERN-SPS in a reasonable amount of time. 

Finally, we stress that the present study is essentially a theoretical one. All the 
"real" work has still to be accomplished in order for one of these options to become 
reality. 



20 



References 



[1] Y. Fukuda et al. (Kamiokande Collaboration), Phys. Lett. B 335, 237 (1994). 

[2] Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81, 1562 

(1998) [|arXiv:hep-ex/9807003|| . 

[3] B. Pontecorvo, J. Expt. Theor. Phys. 33, 549 (1957) [Sov. Phys. JETP 6, 429 
(1958)]; B. Pontecorvo, J. Expt. Theor. Phys. 34, 247 (1958) [Sov. Phys. JETP 
7, 172 (1958)]; Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 
(1962) 870; B. Pontecorvo, J. Expt. Theor. Phys 53 (1967) 1717; V. Gribov and 
B. Pontecorvo, Phys. Lett. B 28, 493 (1969). 

[4] "E362 Proposal for a long baseline neutrino oscillation experiment, using KEK-PS 
and Super-Kamiokande", February 1995. 

H. W. Sobel, proceedings of Eighth International Workshop on Neutrino Tele- 
scopes, Venice 1999, vol 1 pg 351. 

[5] E. Abies et al. [MINOS Collaboration], "P-875: A Long baseline neutrino oscilla- 
tion experiment at Fermilab," FERMILAB-PROPOSAL-P-875. 
The MINOS detectors Technical Design Report, NuMI-L-337, October 1998. 

[6] ICARUS Collab., "ICARUS-II. A Second-Generation Proton Decay Experiment 
and Neutrino Observatory at the Gran Sasso Laboratory", Proposal Vol. I & II, 
LNGS-94/99, 1994. ICARUS Collab. "A first 600 ton ICARUS Detector Installed at 
the Gran Sasso Laboratory", Addendum to proposal, LNGS-95/10 (1995). F. Ar- 
neodo et al. [ICARUS Collaboration], "Cloning of T600 modules to reach the design 
sensitive mass", LNGS-EXP 13/89 add. 2/01, Nov. 2001 and references therein. 
Available at fittp://www. cern. ch/icarus) . 

[7] K. Kodama et al. [OPERA Collaboration], "OPERA: a long baseline v T appearance 
experiment in the CNGS beam from CERN to Gran Sasso", CERN/SPSC 99-20 
SPSC/M635 LNGS-LOI 19/99. 

[8] G. Acquistapace et al., "The CERN neutrino beam to Gran Sasso (NGS)", Con- 
ceptual Technical Design, CERN 98-02 and INFN/AE-98/05 (1998). 

[9] R. Baldy, et al. "The CERN neutrino beam to Gran Sasso (NGS)", Addendum to 
report CERN 98-02, INFN/AE-98/05, CERN SL-99-034 DI and INFN/AE-99/05 

(1999) . 

[10] M. Apollonio et al. (CHOOZ collaboration), Phys. Lett., B466 (1999) ,415 



[11] Y. Itow et al., "The JHF-Kamioka neutrino project", KEK report 2001-4, [hep- 
ex/01060T9| . 



21 



[12] A. Rubbia and P. Sala, "A low-energy optimization of the CERN-NGS neutrino 
beam for a theta(13) driven neutrino oscillation search," |arXiv:hep-ph/0207084| . 



[13] A. Fasso, A. Ferrari, J. Ranft and P.R. Sala, "FLUKA: Status and Perspectives 
for Hadronic Applications" , Proc. of the Monte Carlo 2000 Conference, Lisbon, 
23-26 October 2000, A. Kling et al. eds. , (Springer, Berlin, 2001) 

[14] A. Fasso, A. Ferrari, J. Ranft and P.R. Sala, "Electron-photon transpor t in 
FLUKA: status", Proc. of the Monte Carlo 2000 Conference, Lisbon, 23-26 Octo- 
ber 2000, A. Kling et al. eds. , (Springer, Berlin, 2001). 

[15] D. Harris, "Comparison of Different Detectors with same Beam", 
presented at NNN02, CERN 16-18 January 2002. Available at 
\http : / '/ 'www . cern. ch/muonstoragerings/Nu WorkshopOQ . 

[16] B. Autin et al. (the SPL study group), " Conceptual design of the SPL, a high- 
power superconducting H~ Linac at CERN", CERN 2000-012 (2000). 

[17] M. Shiozawa, XXth International Conference on Neutrino Physics and Astro- 
physics, Neutrino 2002, May 2002 (Germany). 

[18] M. Reiser, "Theory and Design of Charged Particle Beams", (John Wiley Sz Sons 
Inc, New York, 1994). 

[19] J.D. Jackson, "Classical Electrodynamics", third edition, (John Wiley & Sons 
Inc., New York, 1998). 

[20] A. Bueno, M. Campanelli and A. Rubbia, Nucl. Phys. B 589, 577 (2000) 
|aiXiv:hep-ph/0005007|| . 



22 



> 

CD 

O 



CO 




rv Momentum ( GeV/c 



Figure 1: Normalized pion production rate Y^/Ep for various incident proton energies 



E p as estimated with FLUKA. 
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Figure 2: Rescaled energy spectrum of charged current events N v cc{E v )l 'E p for var- 
ious incident proton energies E p (arbitrarly normalized to a baseline L = 732 km). 



24 




Figure 3: Doubly rescaled energy spectrum of charged current events 
N v ^cc{E v ) I 'Ell 'E p for various incident proton energies E p . 
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4.4 GeV beam on a Hg target 




Figure 4: Particle trajectories in the current loops focusing, from a 4.4 GeV proton 
beam on a Hg target. Dotted tracks: p w < 0.5 GeV, dot-dashed tracks 0.5 < p w < 1 
GeV, continuous tracks 1 < p w < 2 GeV, dashed tracks p n > 2 GeV. The thick curve 
is the on-axis magnetic field value (scale on the right). Loop positions are marked on 
top and bottom. 
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4.4 GeV beam 
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Figure 5: Effect of a multiple current loop system on the transverse momentum distri- 
bution, in the case of a 4.4 GeV beam on a mercury target 
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20 GeV beam 
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Figure 6: Particle trajectories in the horn+reflector focusing, from a 20 GeV proton 
beam on a C target. Dotted tracks: p w < 0.5 GeV, dot-dashed tracks 0.5 < p n < 1 
GeV, continuous tracks 1 < p n < 6 GeV, dashed tracks p n > 6 GeV. 
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Figure 7: Effect of a double horn focusing on the transverse momentum distribution, 
in the case of a 20 GeV beam on a carbon target 
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Figure 8: Integrated beam intensity needed to obtain 5 oscillated u e events in a 2.35 
kton detector for Am 2 = 3 x eV 2 and sin 2 (2# 13 ) = 10 -3 for different beam energies 
and baselines. Full symbols correspond to the real focusing systems described in the 
text, open symbols to energy-independent 20% scaling from perfect focusing. 
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E p =2.2 GeV, real focusing 
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Figure 9: Charged current event rate at 120 km as a function of the neutrino en- 
ergy computed with (real) coil focusing for 2.2 GeV proton energy. The dotted line 
corresponds to the oscillation probability (arb. norm.) for a Am 2 = 3 x 10~ 3 eV 2 . 
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E p =4.4 GeV, real focusing 
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Figure 10: Charged current event rate at 120 km as a function of the neutrino 
energy computed with (real) coil focusing for 4.4 GeV proton energy. The dotted line 
corresponds to the oscillation probability (arb. norm.) for a Am 2 = 3 x 10~ 3 eV 2 . 
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Figure 11: Expected 90% CL. sensitivity to sin 2 2^13 mixing angle at 2.2 GeV proton 
energy for the three possible baselines L = 120, 270 and 465 km (thick lines), compared 
to expected results from JHF-SK|TT[ (thin line). 
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Figure 12: Expected 90% CL. sensitivity to sin 2 2#i3 mixing angle at 4.4 GeV proton 
energy for the three possible baselines L = 120, 270 and 465 km (thick lines), compared 
to expected results from JHF-SK|TT[ (thin line). 
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E p =20 GeV, real focusing 
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Figure 13: Charged current event rate at 730 km as a function of the neutrino 
energy computed with (real) horn focusing for 20 GeV proton energy. The dotted line 
corresponds to the oscillation probability (arb. norm.) for a Am 2 = 3 x 10~ 3 eV 2 . 



35 



1 0"' 



fit 



10 



-2 



Three neutrino mixing 



ICARUS, 2,55 kton, 5 years 



IHF+SK, 22.5 



CERN PS++ 20 Ge\ 
L=750 km 
A) 2 10 22 pots/yr 
2 1 21 pots/yr 




on, 5 years 



10' 



10' 



10~ 1 



sin 2< 



13 



Figure 14: Expected 90% C.L. sensitivity to sin 2 2^13 mixing angle at 20 GeV proton 
energy for the CERN-LNGS baseline (L=730 km) (thick lines) for two yearly integrated 
proton intensities, compared to expected results from CNGS r and L.E. optimized (See 
Ref.|2|) (thin lines) and JHF-SK[|nj (dotted line). 
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Figure 15: Expected ratio of charged K + to 7r + production as a function of the primary 
proton energy, for a lm long, 2 mm radius graphite target. The arrows indicate proton 
energies of 2.2, 4.4, 20, 50 and 400 GeV. 
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Figure 16: Ratio of pion decay over the square root of muon decay as a function of 
decay tunnel length, normalized to the value at one pion decay length. 
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Figure 17: Rescaled ratio of muon neutrino flux over the square root of the electron 
neutrino flux divided by the square root of the proton energy (see text) [y^j t/iQ / \fE v 
as a function of the baseline L. 
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